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A compact star with superconducting quark core, the hadron crust and the mixed phase between 
the two is considered. The quark meson coupling model for hadron matter and the color flavor 
locked quark model for quark matter is used in order to construct the equation of state for the 
compact star. The effect of pairing of quarks in the color flavor locked phase and the mixed phase 
on the mass, radius, and period of the rotating star is studied. 
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I. INTRODUCTION 



Neutron/quark stars are the stellar objects produced as a result of supernova explosions. These objects are highly 
compact as their masses may, generally, be one or two times the solar mass and radii of only about 10 — 15 km. The 
neutron stars are normally treated as zero temperature stellar objects because their temperatures, although high, 
are still very low compared to their characteristic energies of excitation. Their outer crust is normally believed to 
be composed, mainly, of nuclei and electrons. Their interior, however, where density is of the order of five to ten 
t-H ■ times the nuclear saturation density, remains to be properly understood whether the central part of the star is either 
composed of quark matter alone, or of mixed matter, or of paired quark matter is one of the subjects of the present 
work. 

Recently, Alford and Rajagopal pj discussed the difficulty of the 2-flavor superconducting (2SC) phase to achieve 
charge neutrality and therefore they claimed that this phase would not be realised in compact stars. Without solving 
the gap equation, they considered strange quark mass as a parameter in their calculations and concluded that the 
color flavor locked (CFL) phase would be favourable against 2SC. In Ref. Q, Steiner, Reddy and Prakash investigated 
the color superconducting phase of quark matter by utlizing the Nambu-Jona-Lasino (NJL) model suplemented by 
diquark and the t'Hooft six fermion interactions. They found that in the NJL model, a small 2SC window does exist 
at very small baryon densities which eventually would be shut by the hadronic phase. Again, the 2SC phase is less 
likely to occur in compact stars because the charge neutrality condition imposes a strong constraint to the quark 
chemical potential. Further, in a recent work [j| suggests that when the density drops low enough so that the mass 
> ■ of the strange quark can no longer be neglected, there is a continuous phase transition from the CFL phase to a 
new gapless CFL (gCFL) phase, which could lead observable consequences if it occured in the cores of neutron star 
However, it now appears that some of the gluons in the gCFL phase have imaginary Meissner masses indicating 
towards an unknown lower energy phase 0, |(| . 

The observations of the fastest pulsars with high degree of compactness and the better theoretical understandings 
of the matter at high densities like that of a compact star with superconducting quark core in a color flavor locked 
phase are compelling theoreticians to come up with different models for the equation of state (EOS) of neutron stars 
and their interiors. 

In the present paper, we are interested in building the EOS for mixed matter of quark and hadron phases. We 
employ the quark-meson coupling model (QMC) 0, Q including hyperons in order to describe the hadron phase. In 
the QMC model, baryons are described as a system of non-overlapping MIT bags which interact through the effective 
scalar and vector mean fields, very much in the same way as in the Walecka model (QHD) 9j. Many applications and 
extensions of the model have been made in the last few years 01 ■ 

While the QMC model shares many similarities with QHD-type models, it however offers new opportunities for 
studying nuclear matter properties. One of the most attractive aspects of the model is that different phases of hadronic 
matter from very low to very high baryon densities and temperatures can be described within the same underlying 
model. This model describes nuclear matter with nucleons as nonoverlapping MIT bags and the quarks inside them 
couple to scalar and vector mesons. For matter at very high density and/or temperature, one expects that baryons 
and mesons dissolve and the entire system of quarks and gluons becomes confined within a single, big MIT bag. 
Another important aspect of the QMC is that the internal structure of the nucleon is introduced explicitly. It is found 
that the EOS for infinite nuclear matter at zero temperature derived from the QMC model is much softer than the 
one obtained in the Walecka model [9j. Also, the QMC model nucleon effective mass lies in the range 0.7 to 0.8 of 
the free nucleon mass, which agrees with results derived from non-relativistic analysis of scattering of neutrons from 
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lead nuclei [llj and is larger in comparison with Walecka model effective mass. 

For the quarkphase, we consider pairing of quarks described by the color-flavor locked (CFL) phase. Recently 
many authors have discussed the possibility of quark matter in a color-superconducting phase, wherein quarks 
near the Fermi surface are paired and these Cooper pairs condense and break the color gauge symmetry |13|. At 
sufficiently high density the favored phase is called CFL, in which quarks of all three colors and all three flavors are 
allowed to pair. 

We organize the paper as follows: In section 2, we briefly discuss i the construction of the equation of state for 
a compact star. The mixed EOS would be built by enforcing appropriate Gibbs criteria and chemical equilibrium 
conditions. In section 3, we present our numerical results and discuss the properties like the maximum mass, radius 
and the period of the compact star, etc. under the influence of the pairing gap, A and the bag pressure, B. 



II. EQUATION OF STATE 

We now discuss the EOS for a compact star with a minimal mathematical details. The details of the QMC model 
can be found in 0, and the details of CFL model can be found in . 

In QMC model, the nucleon in nuclear medium is assumed to be a static spherical MIT bag in which quarks interact 
through the scalar and vector, a, u and p mesons. And these fields are treated as classical fields in the mean field 
approximation. The quark field, 4> q {x), inside the bag then satisfies the equation of motion: 



ip q (x) = o , 



it, d, s, 



(1) 



where m q is the current quark mass and g q , g% and g q p are the quark-meson coupling constants. After enforcing the 
boundary condition at the bag surface, the transcendental equation for the ground state solution of the quark (in s- 

state) is jo(x q ) = (3 q ji(x q ) , which determines the bag eigenfrequency x q , where j3 q — yj (fl q — Rsm*) / (Tl q + Rsm*), 

with Q q — (x 2 + R 2 m* 2 ) 1 / 2 ; m* = m q — g^ao, is the effective quark mass. The energy of a static bag describing 
baryon B consisting of three ground state quarks can be expressed as 
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where Bb is the bag constant and Zb parameterizes the sum of the center-of-mass motion and the gluonic corrections. 
The medium dependent bag radius, Rb is then obtained through the stability condition for the bag. 

The total energy density, e, and the pressure, p, including the leptons can be obtained from the grand canonical 
potential and they read 
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In the above equations 7 and ks are respectively the spin degeneracy and the Fermi momentum of the baryon species 
B. The hyperon couplings, g„B = x a B g<jN, g^B = %uB g^N, gpB — x p b g P N are not relevant to the ground state 
properties of nuclear matter, but information about them can be obtained from the levels in A hypernuclei [l4j . The 
x (j b , x u b and x p b arc equal to 1 for nucleons and acquire different values in different paramctcrizations for the other 
baryons. Note that the s-quark is unaffected by the a and lo— mesons i.e. g% = = . The lepton Fermi momenta 
are the positive real solutions of (k 2 + m 2 ) 1 / 2 = fi e and (k 2 + m 2 ) : / 2 = [i^ = /i e . The equilibrium composition of the 
star is obtained with the charge neutrality condition at a given total baryonic density p = J2b 7 ^/(fr 71 " 2 ); the baryon 
effective masses are obtained self-consistently in the bag model. For stars in which the strongly interacting particles 
are baryons, the composition is determined by the requirements of charge neutrality and /3-equilibrium conditions 
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under the weak processes B\ — > B2 + I + vi and B2 + I — > B\ + vi. After deleptonization, the charge neutrality 
condition yields gtot = 1b 7 k B /(Qir 2 ) + ~^2i =e ^qikf / (3tt 2 ) — , where q B corresponds to the electric charge 
of baryon species B and qi corresponds to the electric charge of lepton species I. Since the time scale of a star is 
effectively infinite compared to the weak interaction time scale, weak interaction violates strangeness conservation. 
The strangeness quantum number is therefore not conserved in a star and the net strangeness is determined by the 
condition of /3-equilibrium which for baryon B is, then, given by /is = bBp n — qBP-e, where ps is the chemical potential 
of baryon B and bs is its baryon number. Thus the chemical potential of any baryon can be obtained from the two 
independent chemical potentials p n and \i e of neutron and electron respectively. 

We next study the equation of state of a compact star taking into consideration the CFL quark phase. We treat the 
quark matter as a Fermi sea of free quarks with an additional contribution to the pressure arising from the formation 
of CFL condensates. 

The CFL phase can be described with the help of the thermodynamical potential which reads [Tlj : 

ttcFLiVq, Ve) = ^quark(Pq) + ^Gs(A*«)Me) + ^l(Pe), (5) 

where fx q — fJ, n /3 and 

^quark(Pq) = — / p 2 dp{p - fig) + — / p 2 ' dp{yj ' p 2 + 7V 2 - fi q ) ^ + B ' ( 6 ) 

71" JO w JO n 

with m u = rrid set to zero, 



u = 2^,-^1 + ^, (7) 

&>GFs(lJ>n, Ue) is the contribution from the Goldstone bosons arising due to the chiral symmetry breaking in the CFL 
phase llB|!il: 

1 / m 2 \ 2 

n GB { N ,» e ) = --f*vl i--f , (8) 



where 



(21-8 ln2)fi 2 3A 2 

£li(n e ) = — (/x e ) 4 /127r 2 , and the quark number densities are equal, i.e., 

+ 2A% 

Pu=Pd = Ps = ~ ■ (10) 

and in the above expressions A is the gap parameter [Til ]. 

We next consider the scenario of a mixed phase of hadronic and quark matter where the hadron phase is described 
by the QMC model and the quark phase is described by the CFL model. In the mixed phase charge neutrality 
is imposed globally i.e. the quark and hadron phases are not neutral separately but rather, the system prefers to 
rearrange itself so that 

X p? p + (i-xVf p + ^ = o (ii) 

where pQ p and p^ p are the charge densities of quark and hadron phases, respectively. \ is the volume fraction 
occupied by the quark phase, (1 — x) is the volume fraction occupied by the hadron phase and p c is the lepton charge 
density. As usual, the phase boundary of the coexistence region between the hadron and quark phase is determined 
by the Gibbs criteria. The critical neutron and electron chemical potentials and the critical pressure are determined 
by the conditions, /if p = pf P = m, i = n,e, T HP = T« p , P HP (p HP ,T) = pQ p (pQ p 7 T) reflecting the need of 
chemical, thermal and mechanical equilibrium, respectively. The energy density and the total baryon density in the 
mixed phase read: 

e = X e QP + (1 - x)e HP + e l , p = XP QP + (1 - x) P HP ■ (12) 
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III. RESULTS AND DISCUSSION 

We start by fixing the free-space bag properties for the QMC model. In the present calculation, we have used the 
current quark masses, m u = m c i = MeV and m s = 150 MeV and R = 0.6 fm for the bag radius. There are two 
unknowns, Zg and the bag constant Bb- These are obtained as usual by fitting the nucleon mass. M = 939 MeV 
and enforcing the stability condition for the bag. The values obtained for Zb and Bb are given in 

La- Next we fit 

the quark- meson coupling constants g a = g^ = 3g^ and g p = g q p for the nucleon to obtain the correct saturation 
properties of the nuclear matter, Eb = e/ p — Al = —15.75 MeV at p = po = 0.15 fm~ 3 , a sym = 32.5 MeV, K = 257 
MeV and M* = 0.774M. We have taken g% = 5.957, g uN = 8.981 and g pN = 8.651. We take the standard values for 
the meson masses, m CT = 550 MeV, m w = 783 MeV and m p — 770 MeV. The meson-hyperon coupling constants used 
in our calculations are x a b — %u b — x p b — v/2/3 and are obtained based on quark counting arguments |l8j | . 

We first plot the pressure as a function of energy density in Figure 1 for different values of gap parameter, A, 
keeping the bag pressure, B 1 / 4 = 210 MeV constant. We see that for lower the pairing gap the higher would be the 
energy density and the pressure of both the transition points between the hadron crust and quark core. And the EOS 
of mixed phase is wider for lower A. Also the EOS for lower A in the mixed region is stiffer but it becomes softer in 
the CFL region. 

We next proceed to calculate the properties of the compact star using the above EOS. The star is assumed to 
be rapidly rotating, relativistic and compact. The details of the model are given by Kamatsu et. al 19] [known as 
Komatsu, Eriguchi-Hachisu (KEH) method] and Cook et. al [2(j ■ The star is assumed to be stationarily rotating and 
hence have axially, equatorially symmetric structures. The metric in spherical coordinates (t, r, 9, </>) can be written 

as, 

ds 2 = -e 2v dt 2 + e 2a (dr 2 + r 2 dd 2 ) + e 2/3 r 2 sin 2 0(d<p - ujdtf (13) 

where the metric potentials, a, (3, v and lo are functions of r and 6 only (geometrized units, c = G = 1). The matter 
is assumed to be a perfect fluid so that the energy-momentum tensor T ab is given by 

T ab = (e+p)U a U b +pg ab , (14) 

where e, p, U, and g ab are the energy density, pressure, four velocity and metric tensor respectively. It is further 
assumed that the four velocity U a is simply a linear combination of time and angular Killing vectors. The details of 
the calculations for solving the field equations are given in Ref. [lflj . 

The radius of the maximum mass star is sensitive to the low density EOS. We have used Baym, Pethick and 
Sutherland [2l| values for pressure and energy at low baryonic densities. In Figure 2, we have shown the mass - 
radius relation for the hybrid star for different values of A in the range — 150 MeV while keeping the bag pressure 
constant, B 1 / 4 = 210 MeV fixed. The maximum mass, M max , is higher for lower A and it shifts towards the low 
radius side as A decreases. The M max for A = 150 MeV is 1.61 Mq and the corresponding radius, R m ax is 18.54 km, 
where as, for A — MeV M max and R ma x are 2.35 Mq and 16.47 km respectively. For an unpaired quark matter 
star [22T |. for the same bag pressure of 210 MeV, we obtained maximum mass and radius, are 2.36 Mq and 16.12 km 
respectively. To compare with, van Kerkwijk et al. have obtained the mass of 1.86 ± 0.32 for an X-ray pulsar Vela 
X-l |23J. The compactness of a star, which is the ratio of M max to R m ax, therefore, decreases with increasing A. We 
observe from Figure 3 that, for A = MeV the ratio turns out to be 0.142 M Q /km, where as, for A = 150 MeV it is 
0.086 M /km. The general relativistic limit for the ratio is, M/R < 4/9 [24[ (i.e. for a uniform density star with the 
causal equation of state P — e) which is much above our range. The lower limit for the ratio found from observations 
of [2!| is 0.115Mo/km. Hence the hybrid star prefers to have low pairing gap to be compact. 

We observe that, as the bag pressure increases, keeping the pairing gap fixed, both the central energy density, e c 
and M max increases but the radius, R ma x of a star decreases. This is in consistent with [2(|. For a fixed A of 100 
MeV, e c = 7.95 x 10 14 gcuT 3 , M max = 1.32 Mq, and R max = 18.8 km for B 1 / 4 = 190 MeV, where as, for B 1 ' 4 = 210 
MeV, e c = 1.29 x 10 15 gem -3 , M max = 2.17M , and R m ax = 17.4 km. But when the bag pressure is kept constant, 
the decrease in A increases e c . Our results thus reflect that the dense cores of the compact stars favor low pairing 
gap and high bag pressures. 

We have also observed that, as A increases the M max decreases unlike the plots in [27j where they have shown an 
opposite trend i.e, the increase in maximum mass with the pairing gap for pure CFL stars with a broader A range 
they have considered. We would like to mention that, our results are not contradicting theirs inasmuch as we also 
have checked for the pure CFL phase and have seen the same trend. From our results, we observe that the stability 
of the hybrid star will not favor the onset of pure CFL phase at its core, with the EOS we have used as the M max 
lies in the mixed region as in 0] f° r ai l the values of A in the range from — 150 MeV for a fixed bag pressure of 
210 MeV. Thus, the CFL effect is reflected on to the stability of the hybrid star only through the mixed phase. 
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In Figure 4, we have plotted the angular velocity versus the mass of the rotating compact star. For a fixed value 
of A and the bag pressure of B 1 ! 4 = 210 MeV, the angular velocity increases with the mass. And we see that as 
A decreases both the angular velocity and the M max increases. For A = 150 MeV, the angular velocity at M max is 
5.5 x 10 3 s -1 , where as, for A = MeV it is 7.7 x 10 3 s -1 . For an unpaired quark matter star also the angular velocity 
obtained is 7.7 x 10 3 s -1 . The angular velocities for the fastest observed pulsars till date, PSR B1937 + 21 and PSR 
B1957 + 20, are 4033 s _1 and 3908s -1 respectively [2jj. We have also seen that the angular velocity increases with a 
decrease in the radius of the compact star. The effect of the paring gap on the period, P = 2%/fl, of the rotating star, 
for a fixed bag pressure of B 1 / 4 = 210 MeV, is shown in Figure 5. As the pairing gap increases the period increases. 
The limit imposed by general relativity on the period of a relativistic compact star should be, P > 0.24 ms [24|. The 
period for the said fastest pulsars till date is 1.55 ms and 1.6 ms respectively [28|,|2!|. The period in our case lies in 
the range of 0.8 to 1.2 milliseconds for the A range of — 150 MeV. The periods predicted for strange stars are in 
submillisecond range j24|. Thus, the stars with low pairing energies are more compact and are the fast rotors. 

Let us now summarize our results. We have studied the compact stars with a mixed matter of hadrons and 
deconfined quarks using QMC model for hadron phase and CFL quark model for quark phase. We conclude here that 
the dense hybrid stars indeed favor low pairing energy and high bag pressures. And these stars are more massive, more 
compact, and are fast rotors with the periods of sub-milliseconds. Also our results clearly show that the hybrid stars, 
of the kind we have considered, possess a quark matter core only as a mixed phase of deconfined quark matter in a 
CFL phase and hadrons. And no pure CFL quark matter core can be possible on the stability grounds. These findings 
have to be deduced more carefully both from the astrophysical observations of compact stars and from laboratory 
measurements of high density matter to confirm or discard the ranges of the said parameters we have considered. 
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FIG. 1: Equation of state for different values of pairing gap and for a fixed bag pressure, B 1/A = 210 MeV. The stiffer(softer) 
EOS in mixed phase becomes softer(stiffer) in CFL phase. 



8 




R(km) 



FIG. 2: The mass - radius relation for different values of pairing gap and for a fixed bag pressure, B 1/4 = 210 MeV. As the 
pairing gap decreases the M max increases and it shifts towards lower radius side. 
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FIG. 3: The degree of compactness, which is the ratio of maximum mass to the corresponding radius of the compact star, 
decreases with an increasing pairing gap while keeping bag pressure, B 1 ^ = 210 MeV, fixed. 



10 




FIG. 4: Angular velocity versus the mass of the rotating compact star for different values of pairing gap and for a fixed bag 
pressure, B 1/4 = 210 MeV. 
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FIG. 5: The period of the rotating compact star as a function of the pairing gap for a fixed bag pressure, B 1//4 = 210 MeV. 
The stars with low pairing gap seem to have sub-millisecond periods. 



